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Abstract 

There is an extensive literature on magnetic gradient induced spin relaxation. Cates, Schaefer and Happer (CSH, [5]), in a 
seminal paper, have solved the problem in the regime where diffusion theory (the Torrey equation [6]) is applicable using an 
expansion of the density matrix in diffusion equation eigenfunctions and angular momentum tensors. McGregor [TO] has solved 
the problem in the same regime using a slightly more general formulation using Redfield theory formulated in terms of the 
auto-correlation function of the fluctuating field seen by the spins and calculating the correlation functions using the diffusion 
theory Green's function. The results of both calculations were shown to agree for a special case, [10] . In the present work we 
show that the eigenfunction expansion of the Torrey equation yields the expansion of the Green's function for the diffusion 
equation thus showing the identity of this approach with that of Redfield theory. The general solution can also be obtained 
directly from the Torrey equation for the density matrix. Thus the physical content of the Redfield and Torrey approaches are 
identical. We then introduce a more general expression for the position autocorrelation function of particles moving in a closed 
cell, extending the range of applicability of the theory. 
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I. INTRODUCTION 

The problem of relaxation in nuclear magnetic resonance due to field gradients has been discussed by many authors 
but continues to be a topic of current research. Recently attention has been focussed on this subject in connection 
with searches for new P,T violating forces mediated by the hitherto unobserved Axion [T], [2]. 

We give a short, very incomplete, summary of how the field developed until now. In 1950 Hahn [3] used his just 
invented spin echo technique to study the effect of translational diffusion on relaxation in nmr. Torrey, in 1953 [1] 
gave a derivation of the effect of translational diffusion that had been alluded to by Hahn. In 1954 Carr and Purcell [5] 
presented a more elaborate method for measuring diffusion constants using relaxation due to translational diffusion 
in an inhomogeneous field with known gradient. Then in 1956 Torrey jB] introduced a specific partial differential 
equation (Torrey equation) describing the effects of diffusion on relaxation. He showed that under conditions when 
diffusion theory was valid the physics was described by adding a diffusion term to the usual Bloch equations. These 
treatments of diffusion did not take into account the effect of the boundaries of the measurement cell. Ten years later 
(1966) Baldwin Robertson [7] gave an approximate solution of the Torrey equation in a relatively small region, defined 
by 2 parallel planes, where the influence of the boundaries was important. Using the method of phase accumulation 
and assuming the phase distribution to be Gaussian, Neuman, in 1973 [8] gave an approximate solution for planar, 
cylindrical and spherical geometries and showed this was in agreement with Robertson's results. 

In 1987, Cates, Schaefer and Happer (CSH) [9J calculated the relaxation for parameters where the diffusion theory 
is appropriate using second order perturbation theory and an expansion in eigenfunctions of the Torrey equation 
applied to the density matrix. At high densities the perturbation theory breaks down and at low densities (A > R) 
the diffusion theory is invalid. (The present work shows how to go beyond this latter limit.) The authors start with 
the equation of motion for the density matrix in the presence of diffusion, (Torrey equation) [B] 

g-iUM + iW, a) 

They consider the deviations of the magnetic field from the volume averaged field as a perturbation and the volume 
averaged field (taken along z) as the unperturbed system. Then, expanding p in the 'eigenpolarizations' of the 
unperturbed problem and carrying out a perturbation expansion in the field variation, taken to be varying linearly 
with position, they obtain a solution valid to second order in the perturbation. 

McGregor [TO] has given a slightly more general treatment based on Redfield's relaxation matrix theory [TT] , as 
presented by Slichter [12] . The starting point of this treatment is the equation of motion for the density matrix 
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expanded to second order in a perturbation ([T2], equation 5.313) 

^ = I (*) , (o)] + [\p* (o) - (*')] , #r (*)] (2) 

where H* represents the deviation of the field from its volume averaged value and starred quantities are expressed in 
the interaction representation with the volume average field considered as the unperturbed system. 

The results show that the relaxation depends on the auto-correlation function of the fluctuating field (frequency 
spectrum of the field fluctautions) as seen by the spins as they move through the measurement cell and the correlation 
function is determined by the diffusion theory Green's function for the case when diffusion theory is valid. 

For high densities, when the boundary conditions do not play a role, the exact solution obtained by Torrey [5] is 
valid. 

Following this work in 1991, Stollcr, Happer and Dyson [13] have shown how to use the exact eigenfunctions of the 
Torrey equation (Airy functions) to get exact solutions in one dimension, de Swiet and Sen [14] have used this and 
other approaches to study a wider range of geometries. Hayden et al (2004) give a nice discussion of the Gaussian 
phase distribution work along with experimental confirmation in a cylindrical geometry |15j . 

McGregor [10] has shown that the results of his Redfield theory treatment are equivalent to those obtained from the 
Torrey equation [9] for the special case of the high pressure limit in a spherical cell. Nevertheless it is illuminating to 
note that the expansion in the diffusion equation eigenfunctions obtained by CSH [9] is in fact the usual eigenfunction 
expansion of the Green's function and hence the results based on the Torrey equation [5J and those of the Redfield 
theory [10] are identical for all cases considered by CSH. We show this in the next section, with details confined to 
an appendix. Thus the physical content of the two approaches are identical in spite of their rather different starting 
points. 

We then show how these results can be applied beyond the diffusion theory limits by giving an analytic expression 
for the trajectory correlation functions valid for a range of pressures wider than that for which diffusion theory is 
applicable. 



II. EQUIVALENCE OF THE TORREY EQUATION AND REDFIELD THEORY RESULTS WHEN DIF- 
FUSION THEORY IS VALID. 



In the appendix we review the calculation of CSH applied to spin 1/2 and using a slightly altered notation. We 
expand the density matrix in the spin 1/2 operators, cro,±i- 
The result for T x , equation (79 1, compare equation (50), CSH: 



T x V 



Re 



(Dkl, - i2n o ) 



(3) 



is seen to contain the Fourier transform of the eigenfunction expansion of the Green' function, equation (83 ) 

3 , {-?>)<t>p> 



0- v 



(Dkj, - iuS) 



(4) 



so that we have (equation 85 eqn.9 in [TO] ) 



L = t 

Ti 2 



dre^ ( [B, (t)] x [B, (t + r)) x + [Br (t)] y [B, (t + r)] y 



(5) 



Similarly the results for T2 (93 1 are also equivalent to McGregor's results (eqn. 10 in [10]) when we take Q in the 
form 



Dkj, 



^ = G (!*•, 1* lo = 0) = J™ drG t[f, if) 



(6) 
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A. Direct solution using Green's function 



As we have shown that the CSH result in terms of diffusion equation eigenfunctions is identical with the McGregor 
result using the Redfield theory and the diffusion theory Green's function it should be possible to derive the result 
starting with the Torrey equation, (equation ( 42 ) in Appendix A) using ( 44 ) 



We expand p as in ( 56 ) 



— = -T Q p + -Tip + DV p 
at i i 



P (-?,t)=Y J °ifi{~?,t) 



(7) 



(8) 



taking the trace with af obtaining 



dt 



i(A«-A 3 )i 



with 



(9) 



(10) 
(11) 



We will treat the sum on the r.h.s. as a perturbation introducing the Green's function for the unperturbed problem, 
Go (x, t) , satisfying 



^-DV'CoM^ (*)*(*) 

and the boundary condition 

it ■ ^G = 

Then we can convert ^ to an integral equation for f[ 

f[{x,t) = /; (0) + [g (x- x',t- t') — Pi (x% 3 f 3 (x',t')e l ^~ A ^'dx'dt' 

j 

which can be solved by iteration being a solution of ^ with the r.h.s. set equal to 0) 

fi(x, t) = jf > + [g (x- x', t~t')—J2 Pi {x% 3 f} 0) (x', t') e^-^'dx'dt' + 

3 

.. [ [g (x- x', t - t') — Y, Pi e^-^'Go (x' - x", t> - t") x ... 

J J %OLi ^ 

Y [T 1 (x'%, jf } (x",t") e l ^- A ^"dx"dt"dx'dt' 



(12) 



(13) 



(14) 



(15) 



If we now operate on this with d/dt and use (12), noting that we will eventually integrate the result over d 3 x so that 
terms containing DW 2 Gq will vanish because of the boundary condition, we find for the second order term: 



ft(x,t) = ±J2l r i( x )k ei{Ai ~ A))t I G (x-x",t-t") 

lOi l / —f J J 

— Fi (x'% k e i(A *- Ak)t "f' k i0) {x",t")dx"dt" 



(16) 
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and averaging over d 3 x, (. 



I jVv (...): 

— f [ d 3 xd 3 x"dt"Y [T 1 {x)] H e^ A '- A ^ t G (x~x",t-t") 
Vioti J J ~i 3 

^[ri (x'% k e^- Ak)t "f' k w (x",t") 



(17) 



To investigate relaxation we set i = k. As we are interested in relaxation of a spatially homogeneous gas we put 



(x", t) — /■ w (t) = const and take it out of the integral since it is the solution of ([9|) with the r.h.s.=0. Then the 
relaxation rate will be given by 



AO) 



= J J d 3 xd 3 x"dt'' Pi (z)]y e^-^Go (a 

a,- L V /Jjfe 



x'',t-t' 



3,k 



J(h 3 -h k )t" 



(18) 



where i = will give 1/T± and i = 1 (+) will give l/Ta. 

Using equations ([77| [7§| [87j [88|) it is easy to see that we obtain eqn.([79) [84} for l/T x and ([90j [9l]) for 1/T 2 . 

Thus direct solution of the Torrey equation ( 42 ) containing a diffusion term, using the conventional second order 
perturbation theory based on the Green's function for the unperturbed equation yields results in agreement with those 
obtained by McGregor [10] by applying second order perturbation theory to the equation of motion for the density 
matrix (Redfield theory), where diffusion theory only enters through the correlation functions of the magnetic field 
and the physical content of the two theories is identical. 



III. BEYOND DIFFUSION THEORY 



Having shown the equivalence of the CSH treatment based on the Torrey equation to the calculation based on 
Redfield theory when diffusion theory is used in evaluating the correlation functions, we widen the range of applicability 
by introducing a form of the correlation function which is also valid when the diffusion theory breaks down, i.e. when 
the condition A c < I no longer holds (A c — vt c is the collision mean free path and L is a typical size of the containing 
vessel) . 



A. Correlation functions for motion in a closed cell. 

Defining a correlation function as 

Rfg(r) = (f(t)g(t + T)) (19) 
with (...) representing an ensemble and time average, we have the following relations |16j 

d 

Rxv (r) = -j^Rxx (t) 
d 2 

Rw(t) = - — R xx (t) (20) 

so the determination of any one will determine the whole family. 

Barabanov et al |17j have calculated the velocity auto-correlation function for particles moving in a closed vessel 
with specularly reflecting walls. The effect of gas collisions are taken into account. The method was initially [T7] 
applied to cylindrical vessels for a case where only the motion normal to the axis is relevant, and then to rectangular 
shaped vessels |18j . The results have been checked by numerical simulations for many cases [17] , [T5] , [15] . The 



function R xv (r) obtained from R vv (r) by means of equation ( 20 ) , has been applied to the study of a false electric 
dipolc moment effect that arises in magnetic resonance experiments in the presence of an electric field |17) , |19) . The 
result can easily be applied to spherical cavities, the only modification being that the distribution of the angle, a, (the 
angle between the trajectory and the normal to the reflecting surface) will be different in the case of a sphere. The 
correlation function, initially obtained for a single velocity, can be averaged over the appropriate velocity distribution. 
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For simplicity we will concentrate on a rectangular vessel in this work. In that case the motions in each of the 3 
directions are independent [18], so we concentrate on one dimension to begin. Equations (27, 36 and 37) of [17] can 
be combined to give (note the r.h.s. of (33) in that paper should be set equal to unity), 



Rw (t) 



8v 



E 



n=l,3,5. 



?/>B ( T ) 



(21) 



where the wall collision time, t w — 2i?sina/u, for particles with velocity v (in the plane of the trajectory), moving 
in a cylinder or sphere of radius R. For the rectangular case we take a = 7r/2 and R — Li/2, (the length of the cell 
along direction Xi) and then t w — Li/vi for particles in a rectangular vessel, moving along direction Xi with velocity 

Vi. 



nirvi 
■qie 



-T)lT 



me 



~V2T 



and 



1 



Vl,2 



vi - m 



(l±«n) 



2r c 

with r Cj the mean time between collisions. We see that u n r w — nn so that 

i>n (t) 



R vv (t) = Sv1 

n=l,3,5.. 



and R vv (0) = vj pi| . 

Using equations (1201) we find 



Rxx{t) 



E 



i 



n=l,3,5.. 

where the constant of integration has been chosen to satisfy R xx (oo) — and we see that (102) 

8L 2 ^ 1 L 2 



#xx(0) 



F - 



7T 4 ^ n 4 
n=l,3,5.. 



12 



= (x 2 ) 



in agreement with McGregor's result ([10], eqn. 24) from diffusion theory. 
If we introduce dimensionless time r' = r/r c and note that 



s„ = \ 1 - 



2n7r 



with /' = Li/\ c where the collision mean free path, A c = ViT c we can write (26) as 



R xx (t') 



L 2 \ 12- 16 



12/ 7T 2 /' 2 



E 



i 



-(1-s„)t'/2 
(1 - S„) 



-(l+ S „)r'/2 
(1 + S„) 



(22) 

(23) 
(24) 

(25) 



(26) 



(27) 



(28) 
(29) 

(30) 



n=l,3,5.. 

Note that s n can be real or complex representing the transition between diffusive and ballistic behavior. 
Figure 1) shows a plot of R xx {t') / R XX (Q) for various values of V . 
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B. Spectrum of the correlation functions. 



We start with the velocity auto-correlation function equation (21) and take the Fourier transform of (22) using the 
definition of Fourier integral used by McGregor [10] : 



so that 



tp n (w) = 2- 



^ n (r) e=^ T dr 



1 



and following (21) 



M 

^2 



E 



ra=l,3.5.. 



Ipn (g) 
W 2 , 



(31) 



(32) 



(33) 



Then the spectrum of the position auto-correlation function, G xx (u), which determines the relaxation is given by, 
following Q 



y>H _ i6«? 

\w) — _ 9 — ~2 — 



16L 4 



XVjTT 6 

n=l,3,5.. 



E 



T «' rc n=l,3 ) 5. 



1 



(C2- W 2^ +W 2 /T 2 



1 




l 






n e 


(-) 2 -l" 

V n7r / 


2 

+ 


uj'l' 


2 






(71TT)' 2 





(34) 
(35) 



where the last equation is written in terms of a normalized frequency, oj' — ujL/vi and length, V = L/X. 

Figure 2) shows 5 (a/, /') = G xx (u) /G xx (0) as a function of w' with V as a parameter. Note this is for a single 
velocity. Averaging over the velocity distribution is straightforward. 

Taking the limit of (35) for u/ << 1/n and reintroducing w, we find 



G xx (w) = 16£>i 



E 



6 



n=l,3,5.. 0J Z + 



(36) 



where D\ = ViX is the diffusion constant for one dimension. This is the Fourier transform of the diffusion theory 
Green's function for this problem as obtained by McGregor [10], eqn. (24). For high frequencies, assuming I' large, 
we obtain (neglecting the 1 in the denominator in (35): 



G xx (ui) 



2vfr c 



U! 2 (1 + W 2 T C 2 ) 



(37) 



which is identical to McGregor's eqn. (13) [10] for the high frequency limit. In obtaining equation (37) we assumed 

luL 

— » 1 



which of course cannot hold for all n. This means we are not properly accounting for the high n modes, which in 
reality would have a contribution of the form ( 36 1 , which is anyway small for large n, and is responsible for the fact 



that ( 37 ) is independent of the size of the vessel. See the discussion under Fig. 3) in [5] 



IV. DISCUSSION 



The approaches of the two calculations are quite different. We have seen that Cates, Schaefer and Happer [9] solved 
the Torrey equation ( 38 ) by assuming an exponential form for the time dependence of p and expanding the decay 



constant and amplitude in a power series in the fluctuating field, treated as a perturbation. McGregor's approach 
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is based on the Redfield treatment of the equation of motion for the density matrix, eqn. (38 1 without the explicit 
introduction of a diffusion term. Recursion is used to get a second order approximation to this equation and the 
second order term is written in terms of the correlation functions of the fluctuating field components as seen by the 
nuclei [T2]. The diffusion theory is then introduced in the calculation of these correlation functions. Lastly we have 
shown that the same results follow from the recursive expansion of the integral equation, derived by use of the Green's 
function, in the manner of the Born expansion. 

Working out the details of the diffusion theory for a spherical cell McGregor showed that his result is equivalent 
to that of in the high pressure limit with Neuman boundary conditions. We have shown that the two approaches 
give identical results whenever eqn. (38 1 and the perturbation theory is valid, thus clearing up any possible confusion 
as to when one or the other of the two quite different approaches is valid. The physical content of both theories is 
identical. 

We have also presented a more general form of the position auto-correlation function for the case of a rectangular 
cell which is valid beyond the region of validity of diffusion theory. 
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VII. APPENDIX A 

A. Perturbation theory of Cates, Schaefer and Happer [9]. 

The authors start with the equation of motion for the density operator, p, in diffusion approximation: 

% = k [H > P] + m2p (38) 
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The Hamiltonian H is broken up into an main term, and a perturbation : 

H=HW+HW 
H (0) = Ml a z 
H {1) = rjhlti ■ ~$ 



(39) 
(40) 
(41) 



where Q is chosen so that the volume average of Tti is zero and rj is an expansion parameter. (Note: normally 
fl = jB D so here Q 0j i are 1/2 the usual values, Oj i = jB 0! i/2) 
We rewrite [381 as 



~ = i [ft cr z , p] + - ft i • O*', 

at % i L 



DV 2 /9 



(42) 



We will approach the problem using time independent perturbation theory, that is we substitute p = p' e 7 * and 
obtain (p' j=f{t)) 

0= ^ 7 +ir o + DV 2 + |r 1 ^p' (43) 

where r o i are linear operators (we now drop the prime on p' , using p to indicate the time independent solution) 

T p = [n o a z ,p] 



Tip 



We will then expand p in the spherical components of the spin 1/2 operators a 

V± = {&x ± iCy) /2 = <J\.2 
^0 = ® z 



(44) 



(45) 
(46) 



The Oi are seen to have the following properties 

(of*. 



with Mi j2 = ±1, M = 0, and 



Thus 



Tr(afaj) = S t] + 5 jo 5 io = SijCXj 



1 (.7 = 1,2) 
2 (J = 0) 



(47) 
(48) 



(49) 



T ai — 2VL MiCJi 



(50) 



We now follow CSH, 9j, by introducing a perturbation expansion for p(r) and 7 into (43) 



7 = 7 (°) + + ^ 7 < 2) 
As this must hold for any value of r\ we collect terms in equal powers of rj : 

0= Uv + ±r + DvApM 
= (V 0) + Jr + dvA P u + Qr x + 7 «) P o 
= (V 0) + ]r + dv 2 ) P (2) + (jr x + 7 (1) ) + 7 (2 V (0) 



(51) 
(52) 



(53) 
(54) 
(55) 
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We look for a solution in the form 



Substituting into (53) and applying Tr(af-) to the resultant equation yields 



2n o M l + DV 2 ) /< (0) (1») = 



(56) 



(57) 



fi(r) has to satisfy boundary conditions on the surface of the measurement cell. CSH, [§], have taken the von Neuman 
conditions (zero current at the walls) but as they point out the method can be applied to the case where depolarization 
takes place at the walls. In any case equation (57), along with the boundary conditions form an eigenvalue problem. 
The solutions are given by the solution to 



{v 2 + k 2 a )4> a -0 



where the eigenvalues k a are determined by the boundary conditions. Then (57) implies 

7 g> = i2MiQ + Dkl 



(58) 



(59) 



In order to solve for the higher order correction terms to the solution it is usef ul to expand the corrections to /j( r ) 
in a series of the zero order functions, (~7^) = 4>ia, (the eigenfunctons of (58)). 



30' 



(60) 



which form a complete set of functions satisfying the boundary conditions, n — 1 or 2 indicates the order of the 
correction. Thus (54) becomes 



+ 7i,a J ViQa 



= E + 7 2M ^° + ^ V2 ) °* E ^'«£>U + ( 7 r i 



P' 



where we used (50 ). Taking Tr(aJ-) of this last equation yields 



where 



o i • i?, crj 

Making use of the orthogonality of the Q , and taking them to be normalized 

d 3 X(j}*p(/) a = 5j3a 



(61) 



(62) 



(63) 



we multiply (61 ) by 4>* g and integrate over the volume: 



= a g ( 7j (0 2 - i2M g il - Dkj) a% a + \ [Y^ \a) + ^% t 5 aP a g 
= a g ( 7 $ - 7 $) a g % a + \ (fi\ [rj^ |a) + ^M^a, 



using (59), where 



'V 



» J 



(64) 



(65) 
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We note that k a= o = 0, corresponding to a uniform distribution in the cell, is a valid solution and we will seek the 



decay parameters for this mode. Thus we put a = = 0, and i — g in (64) obtaining 

o = J(o|[r 1 ] M |o) + 7 ( 1 ) a J 

= 7$ 



(i) i (P\ [ri] fli< |Q) 

W a f7 (0) -7 (0) 



(66) 
(67) 

(68) 



The matrix element in ( 66 1 is seen to be zero for perturbing fields with a volume average of zero. 



Now we use ( 55 ) to evaluate the second order corrections 

7 r i + 7i,a ) E (T o'p0' a j0',i a + 7i,aCT#a 



30' 



Again taking Tr (c*-) of this equation 



o = « s ( 7 [°2 + 



M g n + DV 2 ) £ + E 7 ^ 



30' 



.(1) 



0' 



= a, 



(tS " 7$) a% a [T l]g . a$ , iQ + a g7 j>« Q 



+ 7i,a*s»^Q! s 



(69) 



(70) 



(71) 



where the last result comes from multiplying by <jfp and integrating over volume. Now taking a — f3 = 0, i = g, we 
find 



3,0' 



(2) 



3,0' 



f (o) (o) A 



(72) 



Our derivation has followed the method of time independent Rayleigh-Schroedinger perturbation theory. The 'states' 
are characterized by two 'quantum numbers' a spin index i and a spatial index a, which can stand for 3 indices, which 
appear when we solve the diffusion equation in 3 dimensions. 



B. Calculation of relaxation times, relation to McGregor's result 



We begin by evaluating (72) for i = 0. Since a a = <j z this will be equal to 1/T\. We have then (j> a =o = 



{/nFikjo) 



d'r^ (T>) [TiU 



(73) 
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We write 



Tt x • ? = ([fix] 

fti • ?, o-„l = 2 ([Ox] - [£)i]_ o-+) 



[r 1 ]_ )0 = 2[n 1 ] + 
[Ti], = -2[n 1 ]_ 



(74) 
(75) 

(76) 
(77) 
(78) 



Thus 



-Y (2) 
li=0,0 



n </5'l [ri],,o 10) 



WO,0 



Va 
+ cc 



1 4 

— = — Re 



It 



(?) 



Mr* 



[fii (?)]+E 

|8' 



" ^ (?) 



d 3 r'd 3 r 



d 3 r'd 3 r 



Now the Green's function for the diffusion equation can be written (see Morse and Feshbach, [ 20 J chapter 7) 



G (?, i|?', t') = u (t - £ ^ (7) ^ (?) e 



-Ufe|(t-t') 



with u (t) , the unit step function. Then the time Fourier transform is (r = t — t') 

G (?, w) = / dre l " T G (?, t|?', t') 
./o 



fVe^V^ (7) ^ (?) 



E 

/3' 



^ (7) (?)' 

(Dfc§, - iw) 



(79) 
(80) 

(81) 
(82) 

(83) 



Comparing to the sum in ( 79 ) we see that we can write («o = 2) 



4 Re 



= 4- 



V 

.Re 



(ire 



/ 2 s > , , 



4 Re / dre 1 ' 



[fii (?)] + G t|?' ; i') d 3 r'd 3 r 
n x (7)] _ [n x (?)] + G (?, t|?', f) Po (r', t')d 3 r'd 3 r 



(84) 



Where w D = "fB and Vo{ r ' 1 1') = l/V is the uniform density of magnetization. Then the joint probability distribution 
of an atom being at r at time t and being at r ' at time t' is G(?,t|?',f)p (r',i') and we see that (following 
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McGregor's notation) 



I roc 

w = Ao = 4Rc / d ^ oT <[«i (*)]- [^i (* + !■)]+) 

dre**^ ([n x (i)] B [Ox (i + r)] a + [fix (*)]„ [fix (t + r)]„ 
2 f dre^ ([fl^ \B X (t + r)] x + [Bi («)]„ [B, (t + r)] y 

J CO ^ 



This is the result of the Redfield theory given as equation (9) in McGregor. 
To calculate T<i we have to evaluate 7^—1 - 

3,/?' 



7+,o " Tj,p> ) ocj 



The non-zero matrix elements arc 



so that 



7 (2) 
'i=+,0 



1 



Re £ 

0' 



[r 1 ] + + = 2p 1 ] z 

[T 1 ] 0i+ = -2[fi 1 ] + 



H 7 +'° 7 o,W 



/3' 



4Re]T 

/3' 



' ■ (0|[n 1 ] z |/3')(/3'|[O 1 ],|0) 
(0|[O 1 ]_|/3')( ) 3'|[n 1 ] + |0) 

(0| [fij, 08' I [Slxl |0) , (0| [O x ]_ |/3') </3'| [n x ] + |0> 



1 _ 4 Re 

3fe ~ ~V~ 



2 Re 



V 



(Dkl) 



0' 



2 f-Dfc|, - iw< 
>p> (?) 



d 3 r'd 3 r+ 



Dkj, 



(?) 



(Dkj, - iuj Q ) 



d 3 r'd 3 r 



From ( 83 ) we write 



4, , (V) <t>p, (?y 



G(?,?',u> = 0) = '£ l 



Dkj, 



/ drG(?,t\?',lf) 
Jo 



so that the first term in ( 90 1 can be written 
7 2 Re 



d.T 



oo 

2 /-oo 



Si r 



T7 



[Si C^)] a G t\?', t + T ) Po (r') dVd 3 r 



<M[Bi (t)] z [B! (t + r)],) 
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in agreement with the second term in equ (10) of McGregor. 



From ( 79 ) we see that the second term is 1 /2Ti so that equation ( 90 ) is equivalent to 

,2 r°° 



i _ i r 



drdB! (VHB, (t + T)] z ) 



which is equivalent to equation (10) of |10) . 

VIII. APPENDIX B, SPIN RELATIONS AND MATRIX ELEMENTS 

[ci, a z ] = ~2(j 1 [cr 2 , a z ] = 2cr 2 [01,02] = °z 



Note 



-4 



n r ^ = [fii\ + a- + (7+ + [n 1 \ z a z 



[rr], =-2[n!]. 



2([Oi] + (7_-[0 1 ]_^ 



£1 1 • O^, (7 

[ri 



-[fii] + £r* + [n 1 ] Jl 2er + 

-2[nj + [r 1 ] +i+ = 2[n 1 ] a 

[ni]_cr,- [OJ z 2c7_ 

2[n x ]_ pri]_ _ = -2 [nj, 



1 


(2n + 


I) 4 


1 




(2n + 


I) 6 


1 





f Q (2n+l) 2 8 
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FIG. 2: Frequency spectrum of the auto-correlation function of fig. 1 as a function of reduced frequency and I' . 
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